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1.0  INTRODUCTION 

Any  real  time  system  that  requires  information  about 
the  ionosphere  simply  cannot  depend  on  manual  scaling  because 
of  both  speed  and  reliability.  For  such  systems  automatic 
scaling  of  ionograms  is  a  necessity.  To  be  successful  an 
automatic  scaling  algorithm  must  be  capable  of  correctly  scal¬ 
ing  both  "easy"  ionograms  recorded  during  quiet  conditions 
with  relatively  low  interference  level,  and  "difficult"  iono¬ 
grams  recorded  during  ionospheric  disturbances  or  high  inter¬ 
ference  level.  After  nearly  five  years  of  development  the 
ULCAR  ARTIST  system  can  perform  the  above  task  reliably  and 
in  real  time  automatically  reducing  ionograms  into  several 
useful  forms. 

The  ARTIST  comprises  hardware  and  software  compo¬ 
nents  and  is  a  part  of  the  ULCAR  Digisonde  256  system.  The 
ardware  is  contained  in  a  10. 5-inch  high  rack  mounted  chassis 
and  consists  of  a  microcomputer  with  an  8086  microprocessor 
and  an  8087  math  coprocessor,  512K  of  RAM,  and  two  disc  drives. 
The  software  consists  of  a  CP/M  operating  system  with  a  Micro¬ 
soft  Fortran  compiler  and  the  ARTIST  software.  Upon  receiving 
a  raw  ionogram,  ARTIST  software  immediately  begins  to  process 
the  data  to  produce  the  virtual  height  trace  for  the  ionogram, 
h'(f),  the  usual  ionospheric  scaling  parameters,  plus  ampli- 
ude,  spread  F  and  Doppler  information.  From  the  autoscaling 
information  the  electron  density  profile  can  be  generated  by 
the  profile-fitting  method  (Reinisch  and  Huang,  1983).  This 
method  is  ideally  suited  for  autoscaled  results  as  discussed 
in  Reference  1.  The  results  of  ARTIST  are  outputted  to  a 
standard  RS232C  port  for  connection  to  a  local  terminal  or 
remote  computer  at  a  central  site.  All  of  the  above  tasks 
are  performed  in  real  time  before  the  next  ionogram  is  ready 
for  processing. 
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Autoscaling  routines  have  been  developed  by  other 
institutions  in  the  past  but  they  work  only  for  the  specific 
types  of  ionograms  they  were  tailored  to,  perhaps  as  a  result 
of  testing  only  on  a  few  ionograms.  To  demonstrate  the  versa¬ 
tility  and  accuracy  of  ARTIST,  the  autoscaling  routine  has 
been  tested  on  some  8000  ionograms  from  Goose  Bay,  Labrador, 
for  the  months  January,  April,  July,  and  September  of  1980. 

This  data  base  is  representative  of  one  year's  data  covering 
all  seasons  and  all  types  of  ionospheric  conditions.  For  the 
same  period  the  autoscaling  results  were  compared  with  man¬ 
ually  scaled  hourly  ionograms  (approximately  2200).  The  com¬ 
parison  of  the  autoscaled  data  to  the  manual  data  was  the 
topic  of  a  test  and  evaluation  report  from  which  the  quality 
of  the  autoscaled  results  is  apparent.  For  example.  Air  Force 
Global  Weather  Central  requirement  for  high- latitude  stations 
for  foF2  is  foF2  ±  1  MHz  in  801  of  the  time;  the  autoscaling 
meets  this  criterion  in  96%  of  roughly  2200  ionograms.  Com¬ 
parison  of  the  ARTIST  h' (f )  traces  with  the  manually  scaled 
h'(f)  traces  for  corresponding  ionograms  shows  very  good 
agreement . 

The  ARTIST  system  outputs  the  following  information 
for  each  ionogram;  virtual  height  trace  h'(f),  and  the  follow¬ 
ing  parameters:  foF2»  foFl,  fminF,  MUF(3000)F,  M(3000)F, 
range  spread  F,  frequency  spread  F,  h'F,  h'F2,  foE,  foEs, 
fminE,  range  spread  E,  frequency  spread  E,  h'E,  h'Es.  The 
maximum  usable  frequency  (MUF)  for  other  than  3000  km  dis¬ 
tances  can  easily  be  determined.  Also  average  echo  amplitude 
values  are  given  for  each  megahertz,  normalized  to  a  100  km 
reflection  height;  the  motion  of  the  ionosphere  is  analyzed 
using  the  Doppler  transition  frequencies.  The  results  for 
one  ionogram  represents  one  data  file  identified  by  station 
code,  coordinates,  date,  time  and  sounder  operation  param¬ 
eters.  '  e  pro r  a  is  constructed  to  facilitate  the  addition 
of  ULCAR's  '  -ecx.'on  density  profile  program. 
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In  the  following  sections  the  ARTIST  software  is 
described  and  the  mathematical  procedures  implemented  by 
ARTIST  are  discussed. 


2.0  FUNCTIONAL  DESCRIPTION 


2.1  Overlay  (BISA,  0,  0)  Program  ARTIST  -  The  Driver  Program 
for  the  ARTIST  Algorithm 

2.1.1  Subroutine  SUMENG 

SUMENG  sums  the  amplitudes  of  a  trace  segment  ac¬ 
cording  to  tagging  (0  or  X  echo).  Input  is  the  starting  and 
last  frequency  in  the  segment  and  the  type  of  echoes  to  sum, 

0  or  X.  The  output  is  the  number  of  points  in  the  sum  and 
the  sum  of  the  specific  amplitudes  on  the  curve. 


2.1.2  Subroutine  INIT 

INIT  initializes  one  or  two  dimensional  arrays  to 
zero.  Input  is  the  array  to  be  initialized  and  its  dimen¬ 
sions,  output  is  the  null  array. 


2.1.3  Subroutine  UNPAK 

UNPAK  extracts  the  amplitude  and  status  information 
as  a  function  of  frequency.  Two  records  of  data  are  unpacked 
at  a  time  where  each  record  contains  ten  frequency  lines  and 
each  frequency  line  contains  128  amplitudes  and  status  tags 
corresponding  to  virtual  heights.  In  each  operation  the  TT 
and  K  values  of  the  preface  are  tested  to  insure  vertical 
ionogram  data.  The  predicted  foE  is  calculated  and  other 
necessary  information,  such  as  day,  hour,  minute,  etc.,  is 
extracted  from  the  preface. 


2.1.4  Subroutine  TSTTK 

TSTTK  checks  that  the  TT  and  K  preface  values  of 
the  first  two  records  of  the  ionogram  correspond  to  vertical 
ionogram  data. 
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2.1.5  Subroutine  SPREF 


SPREF  unpacks  the  preface  and  stores  it  in  the  array 
IPREF.  Odd  number  records  have  the  preface  stored  in  the 
first  24  elements,  even  number  records  use  positions  1  through 
16  and  25  to  32. 

2.1.6  Subroutine  SDATA 

SDATA  unpacks  the  data  of  one  frequency  line  into 
128  amplitudes  and  status  tags  corresponding  to  the  128  pos¬ 
sible  digitized  virtual  heights. 


2.1.7  Subroutine  STAT1 

STAT1  unpacks  the  status  of  data  taken  before  June 
14,  1980  and  stores  the  information  in  the  array  ISTS. 


2.1.8  Subroutine  STAT2 

STAT2  unpacks  the  status  of  data  taken  after  June 
13,  1980  and  stores  the  information  in  array  ISTS. 


2.1.9  Subroutine  NO IS  (Figure  1) 

The  noise  threshold  for  each  frequency  is  determined 
and  all  amplitudes  below  the  threshold  are  set  to  zero.  Each 
frequency  line  is  divided  into  two  regions,  an  upper  and  lower 
64  range  bins  and  the  noise  threshold  is  found  separately. 

The  noise  threshold  is  determined  by  distributing  the  32  pos¬ 
sible  digitized  amplitudes  into  16  levels  and  counting  the 
occurrence  of  each  level  in  the  region.  This  gives  a  distri¬ 
bution  curve  of  occurrence  vs  level.  The  noise  level  is  de¬ 
fined  as  the  first  level  above  the  most  probable  level  having 
at  least  half  the  number  of  occurrences  as  the  most  probable 
level.  The  lower  noise  level  of  the  two  regions  is  selected 
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Figure  1.  Subroutine  NOIS  Flow  Chart 
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as  the  noise  threshold.  This  level  overwrites  the  amplitude 
of  the  first  range  bin  for  future  reference. 


2.1.10  Subroutine  PRNLN 

PRNLN  prints  out  the  data  one  frequency  at  a  time, 
either  the  amplitude  (ID  =  1)  or  the  polarization  status  (ID 
2),  0,  X,  or  B  (oblique)  may  be  printed.  This  routine  is 
mainly  for  diagnostic  use. 


2.1.11  Subroutine  FLM  (Figure  2) 

FLM  marks  the  leading  edge  of  a  pulse.  For  each 
frequency,  the  position  of  the  first  maximum  amplitude  of  an 
0-echo  between  two  zero  amplitudes  is  defined  as  a  first 
local  maximum  (FLM).  There  can  be  unlimited  FLM’s  in  one 
frequency  and  each  is  flagged  by  a  negative  sign.  Isolated 
0-pixels  are  ignored. 


2.1.12  Subroutine  CHKEC2 

This  routine  eliminates  the  double  echoes  of  E  or  F 
traces.  For  each  frequency  the  double  echo  FLM  is  located 
and  a  window  (five  bins  high  for  E  and  nine  bins  high  for  F) 
is  placed  over  the  FLM.  The  corresponding  first  order  echo 
height  is  evaluated  and  a  similar  window  is  placed  on  it. 

The  amplitudes  in  each  window  are  summed  and  compared.  If 
the  sum  from  the  double  echo  window  is  the  smaller,  all  FLM’s 
in  it  are  deleted  (the  data  is  reset  to  a  positive  value). 

The  comparison  is  only  made  when  the  first  order  echo  window 
has  at  least  one  0-echo. 
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Figure  2.  Subroutine  FLM  Flow  Chart 
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2.1.13  Subroutine  SUMAMP 


The  amplitudes  within  a  specific  range  of  a  particu 
lar  frequency  given  as  input  are  summed  and  the  result  is  re¬ 
turned  to  the  calling  program. 


2.1.14  Subroutine  TCENTR  (Figure  3) 

This  routine  locates  the  center  of  the  F-trace 
where  the  echoes  with  significant  amplitude  are  mostly  clus¬ 
tered,  usually  in  the  flat  portion  of  the  trace.  This  is 
accomplished  by  grouping  three  abutting  ranges  together  and 
summing  the  FLM  amplitudes  over  the  first  100  frequencies. 

The  heights  of  the  three  highest  amplitude  sums  are  recorded. 
A  five  frequency  by  three  range  bin  window  is  then  moved 
along  the  frequency  axis  at  the  three  selected  heights  in 
search  of  the  maximum  FLM  amplitude  sum.  The  three  resulting 
3*5  windows  are  inspected  and  the  one  having  the  greatest 
0-echo  amplitude  sum  is  the  main  trace  center. 


2.1.15  Subroutine  SWAP 

The  input  consists  of  four  variables,  IA,  IB,  IC, 
ID,  the  routine  swaps  the  values  of  the  first  two  with  each 
other  and  the  last  two  with  each  other  (i.e.  IA  IB,  IB  -*■  IA 
IC  *  ID,  ID  ♦  IC). 

2.1.16  Subroutine  PLOTRC 

This  routine  allows  the  plotting  of  ionograms  on  a 
line  printer.  Several  modes  are  available:  raw  ionograms 
with  noise  suppressed,  ionograms  with  one  of  the  following; 
FLM's,  baseline,  trace  with  amplitude,  trace  with  tagging,  or 
trace  with  Doppler  level. 
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Figure  3.  Subroutine  TCENTR  Flow  Chart 
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2.1.17  Subroutine  FILLZ 

FILLZ  is  called  with  three  variables  II,  12,  IT. 
The  routine  sets  the  array  ITRCF(I ,  IT)  equal  to  zero  from 
position  II  to  12. 


2.1.18  Subroutine  PFOE 

PFOE  evaluates  the  predicted  foE  value  for  Goose 
Bay,  Labrador,  for  the  year  1980  using  a  table  in  terms  of 
day  of  the  year,  hour  and  minute  of  the  day. 


2.1.19  Subroutine  SRANGE  (Figure  4) 

The  average  range  spread  of  the  trace,  excluding 
the  last  five  frequencies,  is  calculated.  At  each  frequency 
a  four  range  bin  window  is  placed  over  the  identified  trace 
point  and  the  amplitude  in  the  window  summed.  The  window  is 
moved  up  the  height  axis  in  one  bin  steps  and  the  amplitude 
recorded  at  each  step.  The  amplitude  at  each  step  is  compared 
to  the  trace  point  amplitude  until  a  12  dB  drop  is  observed. 
The  range  spread  is  given  by  the  height  difference  of  the  win¬ 
dow  base  position  from  the  trace  point  range.  The  average 
spread  is  calculated  by  the  following  formula 


SPRED  = 


topF 

l  Ah'(f) 

f =minF _ 

topF-minF+l 


where  topF  =  critical  frequency  -  four  frequency  bins.  It  is 
assumed  that  the  height  increments  in  the  ionogram  are  S  km, 
and  the  frequency  increments  0.1  MHz. 


2.1.20  Subroutine  SFREQ  (Figure  5) 

The  routine  evaluates  the  frequency  spread  of  the 
ionogram  defined  as  the  difference  between  fxl  and  fxF2.  To 
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Figure  4.  Subroutine  SRANGE  Flow  Chert 
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evaluate  fxl  we  start  at  fxF2,  use  the  X-hyperbola  as  the 
range  frame  and  sum  the  amplitude  of  all  echoes  in  the  frame. 
The  hyperbola  is  also  shifted  up  in  25  km  steps  to  locate  the 
maximum  amplitude  sum.  The  maximum  sum  is  compared  to  the  0- 
hyperbola's  sum  to  check  for  a  6  dB  drop  per  pixel.  (For 
severe  spread  conditions  where  the  baseline  average  range 
spread  is  more  than  50  km,  a  12  dB  drop  is  used.)  If  the 
amplitude  criteria  is  fulfilled  we  designate  the  cusp  frequency 
fxl.  If  the  criteria  is  not  met  we  increase  the  cusp  frequency 
of  the  search  hyperbola  by  one  step  and  repeat  the  procedure 
until  the  amplitude  criteria  is  fulfilled. 


2.1.21  Subroutine  DROPAMP 

This  routine  is  called  by  the  SFREQ  and  FDOWN  rou¬ 
tine,  it  calculates  the  magnitude  of  the  trace  amplitude  and 
sets  the  limit  of  the  dB  drop  to  be  used  by  the  routines. 


2.1.22  Subroutine  OUTERR 

For  some  ionograms  we  cannot  determine  foF2  due  to 
bad  data  or  absorption.  The  routine  OUTERR  finds  these  cases 
and  outputs  999  for  the  foF2  frequency  and  the  h’(f)  height 
of  the  ionogram  to  indicate  no  reasonable  trace  could  be 
found . 


2.1.23  Subroutine  SORT 

This  routine  sorts  an  input  array  in  ascending 

order. 
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2.2  Overlay  (BISA,  2,  0?  -  The  F-Region  Overlay  Program 
F-Region 

This  is  the  driver  program  for  constructing  the  F- 
region  trace  and  extracting  F-region  parameters. 

2.2.1  Subroutine  RASE  (Figure  6) 

This  routine  constructs  a  baseline  which  outlines 
the  trace  by  following  its  leading  edge.  The  coarse  baseline 
end  (ITEND)  is  determined  as  the  highest  frequency  with  FLM's 
on  five  successive  frequencies.  The  baseline  is  then  gener¬ 
ated  starting  at  the  trace  center  where  a  5  f  by  35  bin  win¬ 
dow  is  placed  and  the  median  range  of  all  center  window  FLM's 
is  determined.  The  range  difference  between  the  highest  and 
lowest  FLM  is  also  recorded.  These  two  values  are  used  to 
set  the  window  size  for  searching  for  FLM's  for  the  next  fre¬ 
quency  (the  range  of  the  preceding  frequency  is  used  if  no 
FLM  is  found).  This  process  continues  toward  higher  frequen¬ 
cies  until  no  FLM's  are  found  for  five  consecutive  frequencies 
after  ITEND,  where  it  is  assumed  the  end  of  the  trace  is  en¬ 
countered  and  the  baseline  point  is  set  equal  to  the  highest 
bin.  Starting  again  at  the  center,  the  process  is  repeated 
toward  the  lower  frequencies  until  the  predicted  foE  is 
reached  or  no  FLM  is  found  on  the  next  five  frequencies. 


2.2.2  Subroutine  MED 

This  routine  is  called  by  BASE  to  find  the  median 
range  of  the  FLM's  in  the  search  window.  This  is  done  using 
a  simple  sort  routine. 


2.2.3  Subroutine  DF0F2  (Figure  7) 

It  determines  the  critical  frequency  by  hyperbolic 
fitting.  Two  parallel  curves  denoting  0-trace  and  X-trace 
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Try  3  ways  (K)  of  selecting  points 
£  do  fitting  between  trace  center 
and  baseline  end.  (I)  (If  baseline 
range  spread  >  50  km  use  all  points) 


Use  3  baseline  pt : 
PI  =  1TEN0  -  1+1 
P2  =  PI  -  k 
P3  =  P2  -  k 
To  construct  a  hy¬ 
perbola  which  is  of 
equal  distance  from 
the  3  pt 


Construct  a  x-hyperbola 
which  is  separated  from 
above  o-hyperbola  by 
1/2  fH 


Sum  o+x  amplitudes  on  hy¬ 
perbola  belt  starting  from 
the  frea.  where  hyperbola 
increases  in  slope  and  is 
on  or  above  baseline 


X  If 

hyperbola  sunr 
basel ine/ 
n.  sum?/' 


^  If  \ 

hyperbola  sum 
\max  imum^/^ 


Record  FoF2,  FxF2 
£  hyperbolic  curves 


Replace  basel in e| 
by  o-hyperbola  ' 


.fitting  tried 


Return 


Figure  7.  Subroutine  DF0F2  Flow  Chart 


are  constructed  with  the  frequency  difference  equal  to  1/2 
the  gyrofrequency .  The  O-hyperbola  is  formed  by  three  points 
on  the  smoothened  baseline.  The  selection  of  the  points  is 
described  in  section  5.0.  The  best  fitted  traced"  are  those 
which  give  the  maximum  summation  of  amplitude  in  the  curve 
belts.  To  avoid  selection  of  an  incorrect  trace  the  follow¬ 
ing  condition  is  employed.  The  maximum  sum  must  be  greater 
than  the  corresponding  amplitude  sum  for  the  baseline  between 
the  same  starting  frequency  and  the  baseline  end. 

2.2.4  Subrout ine  SUMBS 

This  routine  is  called  by  DF0F2  with  a  specified 
frequency  range  and  it  performs  the  amplitude  summation  of 
the  baseline  within  the  given  frequency  interval. 

2.2.5  Subroutine  SUMHYP 

This  routine  is  very  similar  to  SUMBS  except  that 
it  evaluates  the  amplitude  summation  of  either  the  0  or  X 
hyperbola  within  the  given  frequency  interval. 

2.2.6  Subroutine  FDOWN 

Subroutine  FDOWN  lowers  the  extracted  F-trace  to 
follow  the  leading  edge  of  the  pulse.  The  amplitudes  of  the 
identified  trace  are  summed  according  to  the  window  sizes  de¬ 
fined  by  each  two  adjacent  frequencies  (the  window  range  is 
determined  by  the  trace  point  height  and  the  trace  point 
height  of  the  preceding  frequency).  The  trace  is  lowered  in 
1  bin  steps  and  an  amplitude  summation  is  done  at  each  step 
and  compared  to  the  starting  (maximum)  sum.  This  is  repeated 
until  a  6  dB  per  point  drop  in  amplitude  is  found  or  until 
the  15  km  maximum  drop  is  reached.  The  trace  is  now  set  at 
the  new  value  and  the  size  of  the  drop  (i.e.  number  of  bins) 
recorded. 


2.2.7  Subroutine  SMOOTH 


This  routine  smooths  the  coarse  baseline  by  fitting 
straight  lines  to  segments  of  the  baseline.  Starting  at  the 
trace  center,  three  trace  points  are  used  to  construct  a 
straight  line  that  is  equidistant  from  the  starting  three 
points.  This  line  gives  us  the  starting  three  smooth  points. 
Using  these  three  points  and  three  adjacent  baseline  points 
we  construct  all  straight  lines  using  combinations  of  any 
three  of  our  six  points  seeking  the  line  with  greatest  ampli¬ 
tude.  This  line  is  used  to  adjust  the  height  of  the  first  of 
the  three  coarse  points,  giving  us  another  point  for  the 
smoothened  baseline.  The  six  frequency  interval  is  shifted 
down  by  one  frequency  and  again  we  have  a  frame  with  three 
smoothened  and  three  coarse  points  from  which  we  smooth  the 
point  closest  to  the  three  smoothened  points.  This  is  re¬ 
peated  until  the  beginning  of  the  baseline  is  reached. 


2.2.8  Subroutine  OBLIQ 

In  this  routine  the  maximum  usable  frequency  (M’T 
(3000))  and  M(3000)  factor  for  the  ionogram  are  calculator}. 
The  vertical  ionogram  trace  is  transformed  to  the  oblique 
ionogram  by  multiplying  each  frequency  with  the  transmission 
factor  M(h* ) 

fQb  =  MCh»)  x  fv. 

The  MUF(3000)  value  is  obtained  as  the  highest  frequency  in 
the  oblique  ionogram,  and  the  M(3000)  factor  is  found  by  di¬ 
viding  MUFC3000)  by  foF2. 

2.2.9  Subroutine  SM0TH2  (Figure  8) 

This  is  similar  to  SMOOTH  except  that  smoothing  is 
done  from  the  trace  center  toward  increasing  frequencies  to 
the  baseline  end. 
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2.2.10  Subroutine  MINF 


MINF  determines  the  minimum  frequency  of  the  F- 
trace.  Starting  at  the  frequency  where  the  five  previous  fre¬ 
quencies  had  no  FLM's  a  3  *  3  window  is  moved  along  the  trace 
and  the  amplitude  sum  is  performed  at  each  step.  When  a  win¬ 
dow  is  found  that  has  at  least  40%  of  the  amplitude  of  the 
trace  center,  the  frequency  of  the  first  0-echo  in  that  win¬ 
dow  is  designated  fminF. 


2.2.11  Subroutine  MINHT 

Given  a  starting  frequency  of  the  search  the  mini¬ 
mum  height  of  the  trace  (FI  or  F2)  can  be  found.  For  h'F  the 
search  range  is  between  the  beginning  of  the  trace  to  the 
trace  center,  for  h'F2  the  range  is  from  foFl  to  foF2.  Con¬ 
sidering  three  trace  points  at  a  time,  the  routine  selects 
the  minimum  height  of  the  trace  by  applying  two  criteria;  the 
average  height  of  the  three  points  must  be  a  minimum  and  the 
height  differences  among  the  three  points  must  be  the  minimum. 
The  minimum  height  is  set  equal  to  the  average  of  the  three 
points  meeting  the  above  criteria. 


2.2.12  Subroutine  F0F1F 

Given  the  autoscaled  h’CF)  trace  this  routine  deter¬ 
mines  the  presence  of  an  FI  layer  and  when  present  the  criti- 
kal  frequency  of  the  layer.  The  method  relies  on  the  fact 
that  just  before  the  critical  frequency  the  slope  approaches 
infinity  and  just  after  foFl  the  slope  is  close  to  zero  or 
negative.  The  routine  searches  from  1.4  foE  to  2.0  foE  as 
possible  foFl  points.  At  each  point  we  take  the  four  preced¬ 
ing  points  and  the  next  four  points;  using  linear  regression 
on  the  two  sets  of  five  data  points  (note  the  point  in  ques¬ 
tion  appears  in  both  data  sets)  we  determine  the  slope  before 
and  after  the  point  in  question.  We  apply  a  double  criteria 
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to  the  point,  the  slope  before  must  be  greater  than  50.0  km/ 
MHz  and  the  slope  after  must  be  less  than  20.0  km/MHz.  When 
both  criteria  are  met  by  a  point  simultaneously  the  search  is 
stopped  and  the  point  is  designated  foFl  and  the  presence  of 
the  FI  layer  is  recorded.  If  no  point  in  the  range  meets  the 
criteria  we  indicate  no  FI  layer  is  present. 

2.2.13  Subroutine  SLOPE 

Subroutine  SLOPE  is  called  by  the  F0F1F  routine. 
Given  a  particular  frequency  on  the  trace,  SLOPE  uses  the 
given  point  and  the  four  preceding  points  and  the  given  point 
and  the  four  next  trace  points  to  form  two  sets  of  five  points 
For  each  set  of  points  the  slope  is  evaluated  using  linear 
regression  and  compared  to  a  prechosen  value.  Flags  are  set 
depending  on  the  comparisons  and  the  information  is  returned 
to  F0F1F. 


2.2.14  Subroutine  LINREG 

Given  an  input  array  of  five  frequencies  and  heights 
the  slope  of  the  points  is  obtained  by  linear  regression  and 
returned  to  the  calling  program. 

2.2.15  Subroutine  DEPOLE 

This  resets  all  FLM’s  back  to  positive  values  since 
the  FLM’s  are  no  longer  used  in  the  processing  of  ionograms. 


Overlay  (BISA, 


0)  Program  DOPAM 


This  is  the  driver  routine  for  Doppler  analysis. 
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2.3.1  Subroutine  MEDAMP 

This  routine  sets  up  a  five  frequency  window  at 
each  MHz  on  the  trace  points  for  the  E,  Es,  and  F  regions  and 
finds  the  median  amplitude  of  the  vertical  echoes  occurring 
in  the  window. 

2.3.2  Subroutine  DOPAMP 

This  is  called  from  MEDAMP  or  DOPTRN.  It  determines 
the  median  of  amplitudes  or  Doppler  levels  in  the  defined  win¬ 
dow  according  to  the  variable  JDA. 

2.3.3  Subroutine  CALVEL 

This  routine  calculates  the  ionospheric  velocities 
at  transition  points  according  to  the  formula 

v  a  Doppler  frequency  x  C _ 

vf  "  2  x  Doppler  -transition  frequency 

(C  =  speed  of  light). 

2.3.4  Subroutine  STRAMP 

Called  by  DOPAMP,  this  routine  stores  amplitudes  of 
window  and  normalizes  them  to  the  reflection  height  of  100  km 
before  sorting  to  find  median. 

2.3.5  Subroutine  STRDOP 

Similar  to  STRAMP  except  the  Doppler  levels  are 

stored. 

2.3.6  Subroutine  DOPTRN  (Figure  9) 

This  is  s«.v  up  for  the  E,  Es,  and  F-regions.  The 
routine  finds  transition  points  along  the  traces.  A  three 
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Figure  9.  Subroutine  DOPTRN  Flow  Chart 


frequency  window  is  used  to  smooth  the  recorded  Doppler  levels 
by  setting  the  middle  level  equal  to  the  end  levels  if  one, 
the  end  levels  are  equal  and  two,  the  middle  differs  from  the 
ends  by  one  level,  i.e.  323  becomes  333.  Tests  are  made  for 
transitions  and  if  found  they  are  recorded. 


2.3.7  Subroutine  TRAND 

This  is  called  by  DOPTRN  to  test  for  transitions. 
The  Doppler  transition  is  defined  as  occurring  when  the  Dop¬ 
pler  levels  are  the  same  for  five  or  more  frequencies  before 
and  after  the  transition  point  (e.g.  111111-000000). 


2.3.8  Subroutine  MEDIAN 

This  routine  sorts  an  array  from  STRAMF  or  STRDOP 
and  finds  the  median  value  of  the  array. 


2.3.9  Subroutine  PCHAR 

PCHAR  encodes  the  median  amplitude  for  each  MHz  for 
printing  out. 


2.4  Overlay  (BISA,  1,  0)  Program  EREGON 

This  is  the  driver  program  for  calculating  the  E- 
region  traces  (E  and  Es)  and  parameters. 


2.4.1  Subroutine  ETRACE  (Figure  10) 

This  routine  calculates  the  E-region  trace  and  pa¬ 
rameters  describing  the  E-region,  i.e.  foE,  h’E,  y^,  etc. 

The  method  assumes  the  electron  density  profile  is  adequately 
given  by 
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Figure  10.  Subroutine  ETRACE  Flow  Chart 
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h-h- 

N  =  U  -  (-r~)2}. 
m  yE 

From  this  the  virtual  height  is  given  by 

h(f  ) 

h»(f)  =  /  y’  dh, 

0 

and  neglecting  the  effect  of  the  geomagnetic  field  on  the 
group  index  of  refraction,  p',  we  can  integrate  giving 

h»(f>  =  hE  -  yE  +  7  yE  tn  rpr- 

The  function  h'(f)  is  fitted  to  the  E-region  echoes  by  opti¬ 
mizing  the  three  unknowns  hE  =  height  of  the  peak  of  the  E- 
layer  parabola,  yE  =  half  thickness,  and  fE  =  foE,  i.e.  the 
critical  frequency  of  the  E-layer.  These  three  parameters 
are  determined  such  as  to  maximize  the  average  signal  ampli¬ 
tude  of  the  ordinary  vertical  echoes  traced  out  by  the  h'(f) 
function.  The  search  range  for  fE  is  ±.3  MHz  about  the  pre¬ 
dicted  value  when  the  predicted  value  is  greater  than  2 . 5  MHz 
otherwise  an  extended  search  algorithm  is  used.  The  range  of 
hE  is  from  95  to  180  km  and  for  yE  it  is  from  10  to  40  km. 
Once  the  E  trace  is  found  calls  are  made  to  routines  to  1. 
may  extend  foE  one  frequency  bin,  2.  lower  the  trace  to  lead¬ 
ing  edge  of  pulse,  3.  evaluate  fminE,  and  4.  evaluate  h'ECf). 


2.4.2  Subroutine  PARBOL  (Figure  11) 

Called  from  ETRACE  with  the  parabola  parameters 
FP  =  fE,  IH  =  hE,  IY  =  yE,  PARBOL  evaluates  the  amplitude  of 
the  points  traced  out  by  the  corresponding  h’(f)  functions. 
For  a  given  FP  and  IY  it  creates  17  parabolas  by  considering 
all  IH  in  the  range  95  to  180  km  and  selects  the  one  with 
maximum  amplitude.  When  parabolas  have  the  same  amplitudes 
the  one  chosen  is 
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Figure  11.  Subroutine  PARBOL  Flow  Chart 
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1 


the  one  with  the  lower  FP, 

2.  the  one  with  the  lower  IY, 

3.  the  one  with  the  lower  IH. 

This  routine,  also  described  in  section  5.0,  can  also  be 
called  to  evaluate  only  the  amplitude  of  one  parabola.  Cau¬ 
tion  -  the  routine  is  only  valid  for  5  km  steps  in  the  height 
range  as  is  presently  done  on  the  DGS  128PS  data  we  have 
processed. 


2.4.3  Subroutine  XTENDE 

This  routine  allows  foE  to  be  increased  one  frequency 
unit.  The  reason  this  is  done  is  that  the  maximum  amplitude 
trace  drawn  out  by  the  ETRACE  results  is  often  short .  To  bet¬ 
ter  compare  with  manual  scaling  we  allow  foE  to  increase  one 
frequency  unit  if  the  following  criteria  is  met.  A  five  range 
bin  window  is  placed  (starting  one  bin  lower  than  the  foE 
trace  point)  at  the  next  frequency  (i.e.  foE  +  .1  MHz)  and  if 
at  least  three  of  the  five  echoes  are  0  echoes  foE  is  in¬ 
creased,  otherwise  foE  remains  unchanged. 


2.4.4  Subroutine  ES  (Figure  12) 

This  routine  evaluates  foEs  if  present.  A  3  x  3 
maximum  amplitude  window  in  the  E-region  is  found  by  first 
summing  all  frequencies  to  16  MHz  as  a  function  of  range  giv¬ 
ing  the  maximum  amplitude  three  range  bins.  Next  the  three 
range  bins  are  summed  from  foE  -  1.0  to  16.0  MHz  and  the  max¬ 
imum  3><3  window  is  obtained.  The  3x3  window  is  moved  out¬ 
ward,  past  foE,  one  frequency  step  at  a  time.  At  each  step 
we  sum  the  O-amplitudes  within  the  window  and  apply  a  double 
criteria;  the  3x3  window  must  have  at  least  five  0-echoes 
and  the  amplitude  sum  must  be  greater  than  5  x  (AVAMP-10). 

The  factor  of  five  is  because  of  the  five  0-echo  minimum  and 


29 


Figure  12.  Subroutine  ES  Flow  Chart 
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we  allow  a  10  dB  per  point  drop  from  the  average  amplitude 
(AVAMP)  of  the  E-trace.  The  following  safeguards  appear. 


1.  The  resulting  foEs  must  be  greater  than  foE  +  .2  MHz 
otherwise  foEs  is  set  to  foE. 

2.  The  search  must  move  at  least  two  frequency  bins  before 
we  accept  any  foEs  value. 

3.  The  trace  (Es)  is  allowed  to  rise  by  considering  a  3  *  3 
box  one  bin  higher  at  each  step.  If  the  amplitude  sum 
of  the  upper  box  is  greater  than  the  lower,  we  jump  up 
one  step. 

2.4.5  Subroutine  EDOWN 

This  routine  is  similar  to  FDOWN,  it  lowers  the  E 
trace  to  follow  the  leading  edge  of  the  pulse.  This  is  ac¬ 
complished  by  lowering  the  parabola  (decreasing  hj.)  and  sum¬ 
ming  the  amplitudes  until  the  sum  is  6  dB  per  point  lower 
than  the  maximum  amplitude  trace.  The  lowering  is  stopped 
when  the  amplitude  criteria  is  met,  when  a  20  km  maximum 
lowering  is  met  or  when  the  lowest  portion  of  the  parabola 
(h^-yj.)  reaches  an  80  km  minimum.  The  new  parabola  parameter 
and  the  number  of  range  bins  in  the  drop  are  recorded. 

2.4.6  Subroutine  ECUSP  (Figure  13) 

This  routine  defines  the  search  range  for  foE  when 
the  predicted  foE  value  is  less  than  2.6  MHz.  The  method 
sums  the  range  bins  from  95  to  175  km  as  a  function  of  fre¬ 
quency  and  find  the  maximum  amplitude  three  frequency  wide 
sum.  The  range  is  determined  by  searching  to  the  right  until 
a  three  frequency  wide  sum  is  less  than  .5  of  the  maximum  sum 
and  by  searching  to  the  left  until  a  three  frequency  wide  sum 
is  less  than  .75  times  the  maximum  sum.  The  search  range  is 
returned  to  the  calling  program. 
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Figure  13.  Subroutine  ECUSP  Flow  Chart 
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2.4.7  Subroutine  FEMIN 

This  routine  calculates  fminE  by  summing  all  ampli¬ 
tudes  in  the  height  range  from  95  to  175  km  as  a  function  of 
frequency.  Starting  at  the  first  frequency  bin  three  fre¬ 
quency  wide  sums  are  made  until  one  is  obtained  that  is 
greater  than  100  dB.  The  first  frequency  bin  in  the  sum 
having  an  echo  (i.e.  amplitude  >  100)  is  designated  fminE. 


2.4.8  Subroutine  ESDOWN 

This  routine  lowers  the  Es  trace  in  the  same  fash¬ 
ion  the  F  trace  was  lowered,  i.e.  as  function  of  amplitude. 
For  a  very  weak  Es  trace  no  lowering  is  performed.  The  80  km 
lower  limit  used  when  lowering  the  E  trace  is  used  here  also. 


33 


3 . 0  INPUT 


3 . 1  Punched  Card 


Card 

No. 

Variable 

Name 

Card 

Col. 

Format 

Variable  Description 

1 

Jun  14 

N/A 

Free 

Format 

$  denotes  input  data 
before  June  14,  1980; 

1  after  June  13,  1980. 

2 

Kl,  K2, 

K3 

N/A 

Free 

Format 

End  time  of  input  tape 
Kl  =  Julian  day  of  the 

year 

K2  *  hour  of  the  day 

K3  =  minute. 

3 

NUM 

N/A 

Free 

Format 

Number  of  ionogram  to 

be  started  for 

processing. 

3.2  Magnetic  Tape  or  Permanent  File 

Name  TAPElj  No.  of  Files  1;  Density  556  BPI;  No.  of 
Tracks  7;  Type  Bin.  Number  of  words  in  record  216. 


Txrrrrir* 


4.0  OUTPUT 


The  printed  output  is  intermediate  and  final  re¬ 
sults  of  E,  Es  and  F-layers  in  that  order.  The  intermediate 
values  are  such  things  as:  predicted  foE,  iteration  of  trace 
center  (IC,  JC),  amplitude  sum  at  each  step,  no.  of  bins  the 
trace  is  lowered,  amplitude  drop  in  6  dB  steps,  baseline, 
last  frequency  of  baseline,  F-trace  before  and  after  hyper¬ 
bolic  fitting,  range  and  frequency  spread  of  the  trace.  The 
information  above  is  printed  out  using  free  format  and  fol¬ 
lowed  by  a  formatted  print  statement  to  list  20  preface  char¬ 
acters  and  the  final  results  of  some  of  the  parameters  as 
follows:  (MPREF(I) ,  I  =  1,  20),  FMINF,  F0F2,  HPRMF,  SPREDF, 

ISFF ,  FXI,  FMUF,  FM3000,  SPREDE,  ISEF,  HPRMF 2 .  Format:  2012, 
2F6.1,  2F6.0,  16,  3F6.4,  F6.0,  16,  F5.0.  Finally  Doppler  in¬ 
formation  is  printed  including  Doppler  transition  velocities 
and  frequencies,  normalized  amplitude  (set  to  100  km)  at  each 
MHz  for  each  trace.  The  processing  of  each  ionogram  can  be 
terminated  with  a  plot  (ionogram  printout)  displaying  eval¬ 
uated  E,  Es  and  F-trace,  or  plotting  can  be  suppressed  and 
only  data  listed. 

In  addition  to  the  output  file,  final  results  are 
saved  on  two  tapes.  Tape  2  contains  following  records  for 
each  ionogram. 

1.  Preface  characters  and  parameters,  (MPREF(I),  I  =  1,  32), 
FMINF,  F0F1,  F0F2,  FXI,  FMUF,  FM3000,  FMINE,  FOE,  FOES, 
HOM,  YMM,  HPRMF,  HPRMF 2,  HPRME,  HPRMES,  SPREDF,  SPREDE, 
ISFF,  ISEF,  ID0WNF,  ID0WNE,  IDOWNES.  Format:  3211, 
5F4.1,  F4.2,  3F4.1,  8F4.0,  512. 

2.  Number  of  frequency  lines  and  corresponding  heights  of 
F-trace  in  km  NUMF,  CITRCF(I,2),  I  =  MINFF,  IF0M) . 

Format:  3314. 
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3.  Number  of  frequency  lines  and  corresponding  amplitudes 
of  F- trace  NUMF,  (IHHO(I),  I  =  MINFF,  IFOM).  Format: 
3314. 

4.  Number  of  frequency  lines  and  corresponding  heights  of 
E-trace  in  km  NUME,  <ITRCF(I,3),  I  =  MINFE,  IFOE). 
Format:  3314. 

5.  Number  of  frequency  lines  and  corres ponding  amplitudes 
of  E- trace  NUME,  (IHHX(I),  I  =  MINFE,  IFOE).  Format: 
3314. 

6.  Number  of  elements  in  array  DOP  for  F-trace  KFF. 

Format:  3314. 

7.  Doppler  transition  velocities,  frequencies  and  Doppler 
number  changes  for  F-trace  (DOP(I),  TRAN(I),  IOLD(I), 
NEW(I),  1=1,  KFF).  Format:  10  (F5.0,  F4.1,  211). 

8.  Average  amplitude  of  the  trace  normalized  at  100  km  for 
F-trace  KFF,  (DOP(I),  1=1,  KFF).  Format:  12,  18F5.1. 

Records  6,  7,  8  are  repeated  with  the  same  format  for  E 
and  Es-traces.  Tape  4  contains  only  one  record  for  each 
ionogram,  which  consists  of  32  preface  characters  and 
list  of  parameters  <MPREF(I),  I  =  1,  32),  FMINF,  F0F1, 
F0F2 ,  FXI ,  FMUF,  FM3000,  FMINE,  FOES,  HOM,  YMM,  HPRMF, 
HPRMF2,  HPRME,  HPRMES ,  SPREDF,  SPREDE,  ISFF,  ISEF, 

JDOWNF ,  JDOWNE,  JDOWNES.  Format:  3211,  SF4.1,  F4.2, 
3F4.1,  8F4.0,  512. 
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5.0  MATHEMATICAL  OR  LOGICAL  PROCEDURES 


5.1  Subroutine  FNDFOF2 

The  hyperbolic  function  used  to  fit  the  X-trace 
near  its  critical  frequency  is 

hx  =  hQ  +  — hQ,  a,  b  =  coefficients  (5.1) 

and  the  function  to  describe  the  corresponding  0-trace  is 

hQ  =  h0  +  » »  f*  =  f  +  \  fH.  (5.2) 

y  a+bf 

Three  points  are  used  to  construct  the  hyperbola 
which  is  determined  in  such  a  way  that  the  distances  are  x, 

-x  and  x;  where  x  is  a  positive  or  negative  number.  The 
range  deviations  of  the  three  points  (f^,  h^)  from  the  curve 
are 

hi  “  <h0  +  glET. )  =  i  =  1,  2,  3.  (5.3) 

These  equations  can  be  written  in  the  form 

(h1-hQ-x)a  +  (h^-hQ-x)  f^  b  -  1  =  0 

(h2-hQ+x)a  +  (h2-h0+x)  fj  b  -  1  =  0 

(hg-hQ-x)a  +  (hg-hp-x)  f3  b  -  1  =  0. 

From  (5.4a)  and  (5.4b)  one  obtains 

a  =  1  -  f  *  r  1  .  1  } 

hl“^0“x  K^Hq+x 

^  _  lxl  1 

b  -  -  rpq  (np^x  _  h2-ir0'+x)* 

Elimination  of  a  and  b  leads  to  a  quadratic  equation  for  the 
determination  of  x  (F.  Scheid,  "Numerical  Analysis,"  Schaum's 
Outline  Series,  McGraw-Hill,  p.  289,  1968): 


(5.4a) 

(5.4b) 

(5.4c) 

(5.5a) 

( 5 . bb) 
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At  this  point,  the  X-hyperbola  is  specified  and  so  is  the 
corresponding  O-hyperbola. 

Next  we  are  going  to  discuss  how  the  three  points 
are  selected.  Suppose  there  are  N  points  on  the  baseline 
after  the  trace  center,  the  triple-point  sets  to  which  hyper¬ 
bolas  are  fitted  are  as  follows 


Group  1 


Group  2 


Group  3 


PN-2  PN-1  PN 


PN-4  PN-2  PN 


PN-6  PN-3  PN 


5.2  Subroutine  SM0TH2 


The  function  used  for  baseline  segment  fitting  is 


h  =  a  +  bf. 


(5.12) 


Again  three  points  are  used  to  construct  a  straight  line  such 
that  the  points  are  passed  through  by  it  and  are  equidistant 
from  it.  Equation  (5.12)  becomes 


a  +  bf^  -  h^  =  x 
a  +  bf 2  -  h2  =  -  x 
a  +  bf3  -  hg  s  x. 


(5.13) 


(5.14) 


(5.15) 


By  solving  the  equations,  one  obtains 


(5.16) 


x  s  7  ^h2  “  hl  " 


(5.17) 


/ 


(5.18) 


a  =  h.^  +  x  - 

and  the  straight  line  is  defined. 

If  the  smoothing  is  operating  on  a  point  P(I)  pro¬ 
ceeding  in  the  direction  of  the  higher  frequency,  five  more 
points  P(I-3),  P(I-2),  P(I-l),  P(I+1>  and  P(I+2)  are  also 
used.  All  combinations  of  selecting  three  points  for  fitting 
are  performed  and  the  amplitudes  of  the  six  frequencies  on 
the  line  are  summed.  The  point  on  the  line  yielding  the  maxi¬ 
mum  sum  replaces  the  old  baseline  point. 


5.3  E-Region  Trace  Extraction 

A  method  to  extract  the  O-trace  of  the  E-region  has 
been  developed.  The  procedure  assumes  the  electron  density 
profile  for  the  E-region  can  be  represented  by  a  parabolic 
expression  (see  Davies,  p.  135) 

h-h_ 

N  a  N  U  -  (—£>).  (5.19) 

“  yE 

The  virtual  height  can  be  evaluated  by 

h(f  ) 

hf(f)  *  /  y '  dh  (5.20) 

using  Eq.  (5.19)  gives 

f  +f 

h’(f)  *  hj.  —  yE  ♦  y  y£  in  (jr-y)  (5.21) 

E  E 

where  hE  is  the  height  of  the  peak  of  the  E-layer  parabola, 
yE  is  the  semi -thickness  and  fE  is  the  critical  frequency  of 
the  E-layer,  foE.  In  the  evaluation  of  equation  (5.21)  the 
effect  of  the  geomagnetic  field  on  the  group  index  of  refrac¬ 
tion  y*  is  neglected.  This  assumption  has  no  effect  on  the 
data,  which  is  obtained  experimentally  and  we  conment  that 
Eq.  (5.21)  is  used  as  an  analytical  equation  to  give  the 
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shape  of  the  E-region  trace.  By  varying  the  parameters  we 
get  many  different  parabolic  type  shapes  from  which  we  choose 
the  one  with  maximum  amplitude  per  point.  Thus  while  the 
assumption  of  zero  geomagnetic  field  was  made  in  the  evalua¬ 
tion  of  Eq.  (5.21)  we  feel  it  has  little  bearing  in  the  sub¬ 
sequent  results. 

The  solution  for  the  trace  is  obtained  by  maximiz¬ 
ing  the  average  signal  amplitude  of  the  ordinary  vertical 
echoes  traced  out  by  the  h’(f)  function.  By  varying  the 
three  parameters  hE,  yE,  and  fE  we  search  all  possible  curves 
traced  out  and  select  the  one  with  maximum  amplitude.  The 
computation  involves  evaluating  the  height  of  the  trace  at 
particular  frequencies  given  the  three  parameters.  To  better 
understand  this  procedure  we  make  the  following  substitutions 
let  R  =  f/fE»  hgE  =  hj.  -  yE  =  bottom  of  layer,  and  define 
1 +R 

♦  =  R  in  (}--£■) .  Equation  (5.21)  becomes 

h’(f)  =  hBE  ♦  |  yE  $.  (5.22) 

In  the  actual  execution  of  the  program  h*(f),  h£,  yE,  and 
hence  hgE,  are  all  integers  representing  digitized  positions 
in  the  ionogram,  i.e.  bin  numbers.  To  solve  equation  (5.22) 
to  trace  out  h'(f)  curve  on  the  ionogram  we  are  only  inter- 
sted  in  values  of 

7  yE  ♦ 

that  are  integers,  i.e.  positions  on  the  range  grid.  In  par¬ 
ticular,  given  values  of  hE,  yE,  and  foE  we  wish  to  relate 
frequency  (actually  R  =  f/fp)  to  virtual  height  via  Eq.  (5.22). 
Since  the  virtual  heights  in  Eq.  (5.22)  are  in  terms  of  in¬ 
teger  range  bins  what  we  wish  to  solve  is 

+  e  2-  *  N  (5.23) 

yE 
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where  N  is  h’(f)  -  hQE  (an  integer)  for  all  possible  values 
of  N.  In  the  solution,  values  of  N  between  0.5  and  1.5  are 
designated  as  bin  1,  1.5  to  2.5  as  bin  2,  etc.  so  the  equa¬ 
tion  we  actually  solve  is 


R  tn 


1+R  . 

T=K  - 


(N 


.4999)  *  i_ 

yE 


(5.24) 


for  the  range  of  yE  considered  (1  to  8)  and  for  a  range  of  N 
(1  to  25).  This  includes  semi-thicknesses  from  5  to  40  km 
and  covers  the  region  of  an  ionogram  from  60  to  260  km.  By 
solving  Eq.  (5.24)  we  can  generate  a  table  of  bins  vs  R.  Thus 
when  solving  for  the  E-region  trace  we  take  f  and  fE,  form 
their  ratio  and  immediately  have  h’(f)  from  the  table,  hg, 
and  y£.  To  solve  Eq.  (5.24)  (a  transcendential  eq.)  we  re- 
express  it  in  the  form 


*(R)  =  R  tn  ^  -(M  -  .4999)  * 

ts 


(5.25) 


and  use  the  Newton-Raphson  method  to  find  the  solutions  for 
♦(R)  =  0.  For  this  we  need  ^’(R)  which  is 


d»(R) 


_  1+R  *  R  d  1+R 

s  ln  tx  *  R  a*  tn  rnr 


(5.26) 


using  the  relationship  tn  y  *  {7  3x  we  **ave 
d  1+R  _  1-R  ,  1  .  1+R  , 

m  *n  T=1T  *  TO  (T=K  +  ~~}7} 

which  simplifies  to 


1*1. 

nr +  T=*  - 


1-RS 


(5.27) 


(5.28) 


thus 


tn 


1+R 


2R 


*  '  »  e*<r 

1=5  TIr* 

Newton-Raphson  gives  for  the  n+1  guess  of  R 


(5.29) 
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Vl  =  Rn  " 

For  an  initial  guess  we  consider 


(5.30) 


♦  =  R  in  (5.31) 

Forgetting  the  linear  R  term  we  can  write 

e*  =  ^  -*■  e*  -  e^R  =  1  +  R  (5.31) 

R(1  +  e*)  =  e*  -  1 


or 


^Initial 

Guess 


r~ ■  and  <p  is  given  by  Eq. 
e”+l 


(5.23). 


The  table  of  numbers  appears  in  subroutine  PARBOL, 
we  evaluate  the  ratio  f/fE  and  find  the  position  in  the  table 
that  is  greater  than  the  ratio.  That  position  is  equal  to 
the  number  of  range  bins,  i.e.  h(f)  -  hgE. 


5.4 _ Error  Routines  and  Indications 


IERR 

It  will  be  set 
parameters  are 


is  the  error  flag  used  all  through 
when  the  following  cases  occur  and 
assigned  equal  to  99,9  or  999.  for 


the  program, 
the  output 
indication. 


43 


Internal  and  Recovery  Procedures 

Error  Condition  Checks  External  Indications  or  Action  Required 
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6.0  PROGRAM  RESTRICTIONS  AND  TIMING 


6.1  Restrictions 

1.  The  program  processes  only  the  vertical  ionograms  of 
Digisonde  256  or  128PS  data. 

2.  The  maximum  number  of  frequency  of  an  ionogram  is  160. 

If  more  frequencies  are  preferred,  the  dimension  of  ar¬ 
rays  IA,  ITRCF,  IHHO ,  IHHX,  IHHOM,  IHHXM,  IBSUM,  and  the 
variable  MAXF  have  to  be  changed  to  the  desired  values. 

3.  The  program  requires  13 OK  (octal)  of  central  memory  to 
run. 

4.  The  average  run  time  for  one  ionogram  (14  records)  is 
<2  CPU  second. 
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